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Abstract. We present a Monte-Carlo package for simulation of high-redshift supernova data, SNOC. Optical and 
near-infrared photons from supernovae are ray-traced over cosmological distances from the simulated host galaxy 
to the observer at Earth. The distances to the sources are calculated from user provided cosmological parameters 
in a Friedmann-Lemaitre universe, allowing for arbitrary forms of "dark energy". The code takes into account 
gravitational interactions (lensing) and extinction by dust, both in the host galaxy and in the line-of-sight. The user 
can also choose to include exotic effects like a hypothetical attenuation due to photon-axion oscillations. SNOC 
is primarily useful for estimations of cosmological parameter uncertainties from studies of apparent brightness 
of Type la supernovae vs redshift, with special emphasis on potential systematic effects. It can also be used to 
compute standard cosmological quantities like luminosity distance, lookback time and age of the universe in any 
Friedmann-Lemaitre model with or without quintessence. 
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1. Introduction 

The study of brightness of Type la supernovae (SNe) 
at high redshifts has become one of the most important 
tools in observational cosmology over the last few years 
(Goobar & Perlmutter, 1995; Perlmutter et al., 1999; Riess 
et al., 1998, 2001), giving the first direct observational ev- 
idence for a presently accelerating universe. With the in- 
creased capabilities of observing very high-redshift SNe, 
more comprehensive analysis tools are needed (see, e.g., 
Mortsell et al., 2001b). 

In this work we describe a Monte-Carlo simulation 
package called The SuperNova Observation Calculator 
(SNOC) which produces synthetic samples of SN obser- 
vations that can be used to estimate the accuracy of 
the magnitude-redshift method for measuring cosmolog- 
ical parameters as well as quantifying the possibility of 
constraining extragalactic dust properties (Goobar et ah, 
2002a) or investigating the matter distribution in the 
Universe (Mortsell et ah, 2001a). For example, SNOC can 
be used to quantify the deviations in the Hubble diagram 
of Type la SNe due to gravitational lensing or dust ex- 
tinction along the line-of-sight. Tools are also provided to 
compute the possible contamination of an observational 
sample by core collapse SNe and the likelihood of obser- 
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vation of multiple (lensed) images from individual high-z 
SNe (Goobar et ah, 2002b). 

2. The program 

SNOC is a modular F0RTRAN77 program using a random 
number generator described in Walck (1988). The pro- 
gram produces an ASCII output containing one block of 
data for each synthetic "event" . One event corresponds 
to one SN. Each SN is characterized by its redshift and 
SN type (i.e whether is a Type la, Ibc, Iln, IIP, IIL or 
'SN1987A like'). The luminosity distance corresponding 
to the SN redshift is calculated from the user supplied set 
of cosmological parameters where we assume a standard 
Friedmann-Lemaitre (FL) universe with arbitrary choice 
of "dark energy". In order to estimate effects along the 
line-of-sight, e.g., from gravitational lensing and dust ex- 
tinction, the light beam is traced backwards in time from 
the observer to the host galaxy, taking into account the 
non-homogeneity of the universe and possible intervening 
matter. This procedure is shown schematically in Fig. 1. 
The program sequence is as follows: 

— Initialization 

— Read user-supplied parameters (see next section) 

— Loop: event generator 
1. Select SN type 
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Fig. 1. Schematic view of the ray-tracing technique used in 
SNOC. The light-beam is traced back-wards in time, from 
the observer to the host galaxy taking into account the non- 
homogeneity of the universe and possible intervening matter. 

2. Generate SN redshift, z 

3. Calculate distance and lookback time to z 

4. Calculate unperturbed SN lightcurve in a "red- 
shifted" B-band filter 1 

5. Perform ray-tracing 

(a) Calculate gravitational lensing magnification 

(b) Perform differential dust extinction calculation 

(c) Estimate dimming due to exotics, e.g., photon- 
axion oscillations 

6. Calculate apparent SN brightness in BVPJ J filters 

7. Generate host galaxy properties 

3. The user supplied parameters 

Table 1 in Appendix A summarizes the user supplied pa- 
rameters for the simulation. These include the following 
main items: 

— The number and type of SNe to be generated. The 
SN types supported are la, Ibc, Iln, IIP, IIL and core 
collapse SNe resembling SN 1987a. 

— The redshift distribution to be used. The options in- 
clude: 

1. Uniform distribution within a chosen range 

2. Tabulated redshift distribution 

3. Listed individual redshifts 

4. Gaussian probability distribution function with 
given mean value and standard deviation 

5. Redshift distribution generated from star forma- 
tion rates [Ref. Dahlen & Fransson (1999)] 

— Options for generating the SN lightcurve 

— Cosmological parameters: [Hq, Hm, ^x, ax(z)} 

— Gravitational lensing. The user supplies the rela- 
tive fraction of galaxy halos parametrized as singular 
isothermal spheres (SIS), Navarro- Frenk- White pro- 
files (NFW), compact objects or spheres with uniform 
mass density. 

— Dust extinction parameters: mean free paths and fre- 
quency dependence for host galaxy dust extinction, 
dust in intervening galaxies and in the intergalactic 
medium. 

— Choice of observed filters (K-corrections) 

1 Transmission function satisfying B z (A)=B(y^) 



— Galaxy parameters: luminosity and mass distribution 

— SN search conditions 

4. The event generator 

4.1. Supernova Type and redshift distribution 

The SN type is selected as either individual type (e.g., 'Ia') 
or a group of SNe (e.g., 'core') which generates core col- 
lapse SNe with a differential rate as in Dahlen & Fransson 
(1999) and Dahlen & Goobar (2002). 

Several redshift distributions are provided by the code: 
uniform or Gaussian distributions as well as physically 
motivated redshift distributions according to: 

dN _ Aircd 2 L (z) r SN (z) 

dz ~ H(z)(l + z) 2 1 + z ' U 

where r$N is the rate of SNe in comoving coordinates. 
In SNOC we implemented the models of the evolution of 
the comoving rate of SN explosions in Dahlen & Fransson 
(1999). In these models an additional user supplied pa- 
rameter r indicates the minimum time needed to form the 
binary systems believed to be the progenitors of Type Ia 
SN explosions. Typical values range from 0.3 to 3.0 billion 
years. 

It is also possible to generate SNe with redshifts from a 
user provided list of redshifts or user defined distributions. 

4.2. The luminosity distance and apparent magnitude: 
cosmological models 

The luminosity distance and the lookback time to the 
event are calculated assuming an isotropic and homoge- 
neous universe. In such a FL universe, the scale factor 
of the universe and thereby the distance estimates are 
functions of the Hubble constant Hq and the energy com- 
ponents of the universe (see, e.g., Bergstrom & Goobar, 
1999). 

The energy density components are treated as perfect 
fluids characterized by their equation of state, pi = cupi. 
Besides the mass density fijvf {cum = 0), a second fluid is 
used in the simulation, fix , with a user specified parame- 
ter cxx(z). Thus, SNOC can generate luminosity distances 
for any choice of CIm and fix including the cosmologi- 
cal constant (ax = — 1) or a time dependent equation of 
state parameter as predicted by quintessence. In addition 
to the FL luminosity distance, or filled beam distance, the 
distance in absence of any focusing of the light beam, the 
empty beam distance, is calculated. 

The apparent magnitude m of a SN at redshift z . as- 
suming the cosmology 6 — [Slj^, £lx-> &x(z)]i is given by 

m(9,M,z) = M + 5log w [d' L (d,z)}, (2) 
M = 25 + ill + 51og 10 (c/ff ), (3) 

where M is the absolute magnitude of the SN, and 
d' L = Hodi is the Hq -independent luminosity distance, 
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where H is the Hubble parameter. 2 Hence, the inter- 
cept M. contains the "unwanted" parameters M and 77 
that apply equally to all magnitude measurements [we do 
not consider evolutionary effects M — M(z)\. The H a - 
indcpendent luminosity distance d' L is given by 



d'r = 



I 

H'(z) 



(! + z )7=nT sin (V-^fe I), < 

{i + z)i, n fe = o 

n k > o 



(i + z)^= S mh(vn;i), 

- — &X, 

z dz' 



f(z) = cxp 



Jo 77V)' 
77(z)/77 = 

y/(i + z)3 n M + f(z) n x + (i + z) 2 r» fe , 

l + a x {z')' 



3 / dz 1 ' 
o 



l + z 1 



(4) 

(5) 
(6) 

(7) 
(8) 



where we consider an equation of state with arbitrary z- 
dependence. 

4.3. Ray-tracing 

The ray-tracing is discussed in detail in Sect. 5. Here, we 
summarize its key elements. The light-beam path is di- 
vided into 'cells' corresponding approximately to average 
distances between galaxies. Each cell is characterized by 

1. The mass density profile, the source of gravitational 
lensing 

2. The dust density and differential extinction coefficient, 
treated separately for the actual galaxy in the cell cen- 
ter and the region around it 

3. The average magnetic field and electron density, rele- 
vant in exotic scenarios such as the photon-axion os- 
cillation mechanism discussed in Sect. 5.3 

4.4. Supernova and Galaxy apparent magnitudes: 
K-corrections 

The SN and host galaxy apparent magnitudes through a 
standard (or user specified) broad-band filter is calculated 
using cross-filter K-corrections as described in Kim et al. 
(1996). Thus, the apparent magnitude of a SN observed 
at a redshift z through a filter Y defined through the filter 
function Sy(A) becomes: 



K-correction becomes: 

/ JV(A)g^(A)AdA \ 
K XY - 2.51og 10 ^ v(A)Sy(A)AdA j 

/ j7s N (A)MA)ArfA 
gl ° UF SN (A)e i (A)5y(A(l + z))A ( iA; ' 

where the photometry system is normalized in the Vega 
system, V. 

By default, SNOC calculates K-corrections for K BY , 
where Y includes the standard broad-band filters B, V, R, 
I and J. In order to assess the color corrections from the 
various extinction sources we list the results separately: 

0) K BY : no extinction (e = 1) 

1) K BY : host galaxy extinction 

2) K BY : extinction from galaxies along the line-of-sight 

3) K BY : extinction from intergalactic dust 

4) K BY : differential effects from photon-axion oscillations 

Finally, Kby (no index) shows the combined effect on 
the observed magnitudes from all the considered attenua- 
tion sources. 

Similarly, representative galaxy templates for E, Sbc, 
Scd and Im galaxies from Coleman et al. (1980) arc used 
to calculate the observed colors of the SN host galaxies. 

SNOC generates simulated lightcurve data with a user 
specified time sampling, date of discovery and noise con- 
tamination. For Type la SNc, K-corrections for BVRIJ 
may be calculated starting from up to 19 days prior 
to maximum up to 70 days past maximum in the rest- 
frame system using a spectral template from Nugent 
ct al. (2002). 3 For other types of SNe, approximate K- 
corrections are calculated from a set of (truncated) black- 
body functions, following the approach in Dahlcn & 
Fransson (1999). 

4.5. Simulation of supernova search conditions 

Two samplings of the SN lightcurve are generated ran- 
domly using a user specified time gap. SNOC can thus be 
used to optimize the search strategy for finding SNe at any 
given redshift. The tuning parameters are the broad band 
filter chosen, the limiting magnitude of the search and 
the time between the reference and the discovery images. 
The method can also be used to study the contamination 
from core collapse SNe in the discovery sample (Dahlen & 
Goobar, 2002). 



m Y (z) = M x + 51og 10 [d' L (0,z)} + K XY , (9) 

where the restframe absolute magnitude of the SN is mea- 
sured in the filter X with a filter transmission function 
Sx(X). 

Assuming a spectral template for each SN type, TgN 
and an arbitrary differential attenuation < £j(A) < 1 the 

2 In the expression for M, the units of c and Ho are km s _1 
and km s" 1 Mpc -1 , respectively. 



5. The elements of the ray-tracing section 

5.1. Gravitational lensing 

When simulating observations of distant sources, it is nec- 
essary to take into account the fact that the Universe is in- 
homogeneous. First, the inhomogencitics might influence 
the overall expansion rate of the Universe, the so called av- 
eraging problem. Second, gravitational lensing effects due 

3 The spectral template of day 70 is used for later epochs. 
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to the inhomogcncities will cause a dispersion in the ob- 
served source luminosities. 

In SNOC we employ a method proposed by Holz & 
Wald (1998) for examining lensing effects. In this method, 
the overall expansion is assumed to follow that of the av- 
erage background geometry. That is, in the following we 
will neglect the averaging problem and only consider ef- 
fects from gravitational lensing. 4 

5.1.1. The Method 

The method of Holz & Wald can be summarized as fol- 
lows: First, a FL background geometry is selected, i.e., we 
specify the space-averaged energy density content of the 
Universe, and fix ■ Inhomogeneities are accounted for 
by specifying matter distributions in spherical cells which 
have an average energy density equal to that of the under- 
lying FL model. Thereafter, an infinitesimal beam focused 
at the observer is integrated backwards using the geodesic 
deviation equation by being sent through a series of cells, 
each time with a randomly selected impact parameter with 
respect to the matter distribution in the cell. Finally, the 
resulting area, A, at the source redshift is considered. 

In a homogeneous FL universe, the ray bundle will 
have an area corresponding to the filled beam area, Afl- 
Since the apparent luminosity of a source is proportional 
to the inverse of the area of the beam, we define the mag- 
nification as 



By using Monte-Carlo techniques to trace a large number 
of light rays, statistics for the apparent luminosity of an 
ensemble of sources is obtained. 

As noted in Holz & Wald, some care is needed when 
interpreting results as probability distributions. This is 
due to the fact that individual ray bundles do not cor- 
respond to random source positions. Since the fraction of 
the source sphere being sampled by a beam is proportional 
to the area of the beam, i.e., inversely proportional to the 
magnification, magnified sources will be over-represented 
if not compensated for. This is done by designating a prob- 
ability proportional to the area of the ray bundle to keep 
each event. Further discussions of the statistical weight 
given to random lincs-of-sight can be found in Ehlers & 
Schneider (1986) and Bergstrom et al. (2000). 

In order to validate our implementation of the method, 
we have performed a number of consistency tests. The 
interested reader is referred to Bergstrom et al. (2000). 

5.1.2. Matter distributions 

The matter distribution in the cells can be chosen to be 
in arbitrary fractions of point-masses, uniform spheres, 
singular isothermal spheres (SIS) or the Navarro-Frenk- 
White (NFW) density profile (Navarro et al., 1997). 

4 For general references on gravitational lensing see, e.g., the 
textbook by Schneider et al. (1992). 



Matter distribution parameters such as the scale radius 
of the NFW halo and the cut-off radius of the SIS halo 
are determined from distributions reflecting real galaxy 
properties. 

One of the advantages of the method of Holz & Wald is 
that any mass distribution and number density, including 
possible redshift dependencies, can easily be implemented 
and used as long as the average energy density in each cell 
agrees with the underlying FL model. 

We derive a galaxy mass distribution, dn/dM, by 
combining the Schcchtcr luminosity function (see, e.g., 
Peebles, 1993, Eq. 5.129) with a mass-to-luminosity ratio 

MocL^?. (11) 

For elliptical galaxies in the fundamental plane, we can, 
e.g., choose (3 = 0.25 (see Peebles, 1993, Eq. 3.39). The 
normalization of the resulting galaxy mass distribution, 
i.e., the characteristic galaxy mass, M*, is determined 
from observational constraints on the galaxy number den- 
sity. 

5.1.3. Foreground cluster 

Since all cells have the same average density as the un- 
derlying FL model, the average luminosity will be equal 
to the one obtained in the homogeneous case. To simulate 
the increased gravitational lensing effects due to a specific 
matter distribution along a given line-of-sight, the user can 
specify the properties of a cell at a given redshift to mimic 
a foreground galaxy or cluster. The matter distribution is 
in the form of an isothermal sphere, 

Psis(r) = ~ (12) 

where v is the line-of-sight velocity dispersion of the mass 
particles. Thus, the velocity dispersion sets the concentra- 
tion of the cluster, the mass sets the cut-off radius and the 
radius of the cell sets the size of the image field. 

5.1.4. Multiple images 

Since wc in principle only are able to follow infinitesimal 
light-rays in SNOC, we have to use some approximations 
when trying to get information on multiple image systems. 
The main approximation is that we assume that in cases of 
strong lensing, the effects from one close encounter is dom- 
inant, i.e., the one-lens approximation. With this simplifi- 
cation, we can use the information from the magnification 
to derive quantities for systems with finite separations. In 
order to do this, we need to be able to derive analytical 
relations between different observables, e.g., the magnifi- 
cation and the image separation. Here we show how this 
is done for the case of SIS lenses. 

When running SNOC, we want to answer the following 
questions: 

— Will the particular event correspond to the primary 
image in a multiple lens-system? If this is the case: 
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— What is the flux of the second image? 

— What is the separation between images? 

— What is the time-delay between images? 

Employing the one-lens approximation with a SIS lens, 
we will have multiple imaging whenever Mi > 2 and we 
can relate the magnification of the secondary image, M2, 
to that of the first by 



\fi 2 \ = Mi - 2. 
Furthermore, the image separation is given by 



(13) 



(14) 



where D s and D d are the angular size distances, to the 
source (z s ), the lens {z d ), while D ds represents the angular 
size distance between the lens and the source. The time 
delay for the two images is 



At 



4tt 



(- 

V c 



(1 + z d ) 



Mi - 1 



(15) 



Therefore, in order to compute the quantities of interest, 
we need to pick a lens rcdshift and velocity dispersion 
from some reasonable distributions for every case where 
Mi > 2. The distributions will in general be functions of 
the cosmology, the mass distribution of the lenses and the 
source redshift. The differential probability for multiple 
imaging is in the general case given by 



dft 

dP cx a(z d , z s ) — dM{\ + z d ) 3 dV, 



(16) 



where a(zd,z s ) is the cross-section for multiple imaging 
and n is the comoving number density of lenses, assumed 
to be constant. For SIS lenses, the cross-section is given 

by 



-J D, 



2 

ds' 



(17) 



Since dV cx D^j^dz d and n is independent of lens redshift 



we can use 



dt 



dP{z d ) cx Dj s Dj(l + z d Yj-dz dl 



(18) 



as our probability distribution for z d . The probability dis- 
tribution for v is given by 



, , , dn dM /v\ 4 
dP{v) oc -777-7- - 
dM dv \cJ 



dv, 



dM dv vc/ 

where the velocity dispersion is related to mass by 
- = M^K 

(see Bergstrom et al., 2000, Sect. 7.1). 



(19) 



(20) 



5.2. Extinction by dust along the line-of-sight 

An important effect to consider when observing SNe is the 
absorption of light due to dust. With SNOC, it is possible 
to take the effects from three different extragalactic dust 
components into account: 

1. A homogeneous intergalactic component. 

2. A host galaxy dust component. 

3. Dust in galaxies between the source and the observer. 

Properties such as the exponential extinction scale-length 
and the frequency dependence can be set separately for 
each component. The frequency dependence is computed 
using the parametrization in Cardelli et al. (1989). 

5.2.1. Grey dust 

For a given emission redshift z s , the attenuation Am d at 
observed wavelength A Q due to dust can be written 



-2.5 log 



10 



Am d (z s , A Q ) = 

e -C/ z " p Awt {z)a{\ /(l+z),R v )h(z)dz 



(21) 



where Mdust(-z) is the physical dust density at redshift z, 
a(A, Ry) is the wavelength-dependent attenuation, and 
the cosmology-dependent function h(z) is given by 

h{z) = [H{z){l + z)]-\ 



The reddening parameter Ry is defined by 

A v = R V E(B~V), 
where Ay is the V-band extinction and 

E(B — V) = (B — V) — (B — V)i 



(22) 



(23) 



(24) 



with (B — V)i being the intrinsic (unobscured) color. 

The normalization constant C is related to the overall 
magnitude of the extinction, i.e., the exponential extinc- 
tion scale- length. 5 

In a manner similar to that used in computing gravita- 
tional lensing effects, the integral in Eq. (21) is performed 
numerically by following individual light paths through a 
large number of cells containing galaxies and intergalactic 
dust. Through each cell the background cosmology, the 
wavelength of the photon and the dust density are up- 
dated, and the contributions from each cell added. Note 
that the model is approximately valid also for a patchy 
dust distribution, as long as the scale of inhomogeneitics 
is small enough, i.e., 1/V~N -C 1 where N is the number 
of dust clouds intersected by the light-ray. 

The intergalactic dust density is parametrized by 

Pdust = Pdust( 1 + Z ) i(Z \ where 



/(*) - 



a = const, for z < Zy lr 
(3 = const, for z > zu r 



(25) 



5 This can be fixed, e.g., by demanding that a given cosmol- 
ogy reproduces the observed luminosity distance at z ~ 0.5, 
i.e., that dust extinction can explain the dimming of the 
presently observed SN sample, otherwise attributed to the 
"concordance" cosmology Q,m = 0.3, Qa = 0.7. 
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zii m is a user denned redshift where the intergalactic dust 
distribution is assumed to have changed properties. It is 
also possible to supply a user specified set of values of the 
dust density Pdustfz) and the frequency dependence as a 
function of redshift. 



5.2.2. Host galaxy dust 

The absorption in the host galaxy will be different for dif- 
ferent SN types. This is because core collapse SNe (Type 
II, Ibc etc.) occur in late type, star forming, dusty galax- 
ies whereas Type la can occur in any type of host galaxy. 
In the current version of SNOC, we have used Fig. 6 in 
Matteucci & Recchi (2001) to make a first rough estimate 
of the probability for a Type la SNe (up to z < 4) to occur 
in an early type galaxy as 



p E (z) = Pe(0) + z 



(26) 



where Pe(0) = 0.5 and qi — 0.125. This simplistic ap- 
proach could without doubt be improved on in the future. 

Also, the amount of dust and the SN rate depend on 
the size of the galaxy (i.e., the star content). Using Monte- 
Carlo methods, the galaxy type, dust content, SN position 
and disk inclination are determined and the absorption is 
calculated. 

We model the dust distribution in late type galaxies 
as a double exponential disk 



Pdust — Pdust e 



-r/ro-|C|/Co 



(27) 



where r and C are cylindrical coordinates. We set Co = 0.1 
kpc. In order to determine r , we assume that the amount 
of dust is proportional to the amount of stars in the galaxy. 
If stars are distributed as 



n« = nre 



0,,-r/ro-KI/Co 



(28) 



the luminosity, L, will be proportional to J n*dV oc r 2 , ■ Co- 
if Co is fixed, we have r oc \[L. That is, we set ro = 
where r° = 5 kpc is the scale radius for a L — galaxy. 
We determine p^ ust by setting the exponential scale-length 
at the galaxy center. The disk is truncated at r max = 20 
kpc and Cmax = 3 kpc. 

We let the probability for a SN to occur in a specific 
galaxy be proportional to the galaxy luminosity in order 
to allow for SNe to be more common in star-rich, luminous 
galaxies. Core collapse only occur in the disk, closely fol- 
lowing the dust distribution in Eq. (27), but we do not al- 
low any SNe inside r — 0.3 kpc. We let a fraction (87.5 %) 
of Type la SNe follow the dust distribution in Eq. (27), 
only that we set Co = 0.35 kpc, and the remaining frac- 
tion (12.5 %) occur in a bulge according to the probability 
distribution 



p(R) oc 



1 



R 3 + 0.7 3 



(29) 



where R is the radial coordinate in spherical coordinates 
(expressed in kpc). The bulge is truncated at R max = 3 



kpc. For further reference regarding dust and SNe distri- 
bution in host galaxies, see Hatano et al. (1998). After 
determining the dust distribution and SN position in the 
galaxy we assume a random inclination of the disk and 
integrate the absorption along the path of the light-ray 
through the disk. 

5.2.3. Dust in intervening galaxies 

In each cell, the galaxy type is selected using a Monte- 
Carlo method. We use a parametrization where the frac- 
tion of early type, dust-free galaxies is given by /e(z) = 
/b(0) - <72 * z, where f E (0) = 0.3 and q 2 = 0.05. If the 
galaxy is a spiral, the absorption due to the dusty disk 
is computed (if the light-ray passes it). Equivalent to the 
case of host galaxies, we model the dust distribution in 
intervening spiral galaxies by a double exponential disk 
with random inclination, see Eq. (27). 

5.3. Photon-axion oscillations 

The general structure of SNOC, i.e., the tracing of light 
through successive cells, makes it relatively easy to add 
more refinement to the propagation code. As an example 
we mention here an exotic, hypothetical process that could 
affect light propagation on cosmological scales. 

It has recently been proposed that the observed faint - 
ness of high- redshift SNe could be attributed to the mixing 
of photons with a light axion in an intergalactic magnetic 
field (Csaki et al., 2001). We compute the mixing prob- 
ability using the formalism of density matrices (see, e.g., 
Sakurai, 1995). We define the mixing matrix as 



Aj_ Am cos a 

Q = J Ay Am sin a 

Am cos a Am sin a A m 



(30) 



The different quantities appearing in this matrix are given 

by 



A M 



-3.6 x 10~ 25 
Aj_, 

2 x 10~ 26 



10" 8 cm- 3 . 



cm 



Bo,. 



M a 



10- 9 G/ Vl0 n GeV 



cm 



cm 



(31) 
(32) 

(33) 
(34) 



where -Bo,_l is the strength of the magnetic field perpen- 
dicular to the direction of the photon, M a is the inverse 
coupling between the photon and the axion, n e is the elec- 
tron density, m a is the axion mass and u> is the energy of 
the photon. The angle a is the angle between the (pro- 
jected) magnetic field and the (arbitrary, but fixed) per- 
pendicular polarization vector. The time-evolution of the 
density matrix p is given by 



l5tP = 2^' P ^' 



(35) 
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with initial conditions 



Po = 



| 0' 
\ 




(36) 



Here the three diagonal elements refer to two different po- 
larization intensities and the axion intensity respectively. 
In each cell, we solve the 9 coupled (complex) differential 
equations and update p and Q. This method ensures that 
the full frequency dependence of the effect is taken into 
account, as shown in Mortsell et al. (2002). 

6. Fitting cosmological parameters 

The SNOC package also includes a software tool for cos- 
mology fitting, based on the output files from the Monte- 
Carlo program. In its current state two different types of 
maximum-likelihood analysis can be done using different 
measurement variables. 

Type la supernovae as standard candles The analysis tool 
can be used to fit the cosmological parameters M, Om, 
i ctx ; and parameters related to the quintessence mod- 
els presented in Goliath et al. (2001) and Eriksson & 
Amanullah (2002), by using simulated Type la events 
where the magnitude nii (in the redshiftcd B-band) and 
the redshift Zi are considered to be the measurement vari- 
ables for each event i. The negative log- likelihood function 
will in this case take the following appearance 



C = -Nln 



N 



1 



27TCT 



+ 



~2 XI - 51 °gio d L(^M,^x,a;^) - M} 2 , 



2a 2 



i=i 



(37) 



where N is the number of events, a the statistical uncer- 
tainty for each event, nii and Zi and finally the luminosity 
distance, d' L , is given by Eq. (4). 

It is also possible to use a different, more general, like- 
lihood function that also takes into account the fact that 
the data could be subject to systematic effects from grav- 
itational lensing. This procedure is further described in 
Amanullah et al. (2002). 

Cosmological parameters from lensed supernovae The other 
way to use the analysis tool is to use lensed SNe as is de- 
scribed in detail in Goobar et al. (2002b). In this case core- 
collapse SNe as well as Type la's can be used, and in its 
present state the software allows fitting of the parameters 
H , flM, and ax(= cons tant), by using multiply im- 
aged high-z SNe in spherically symmetric lensing systems 
producing only two or ring- like images. The observables 
are assumed to be the source redshift, the redshift of the 
lensing galaxy the image separation, A0, the time delay, 
At, and the flux-ratio, r/, between the images. 



The negative log-likelihood function for this analysis 
can be expressed as 



r _ sr^ (-Rpxp TV) 2 
^ 2a 2 



(38) 



where W cxp and TZ 1 can be written as 



D d D s 



TV 



D ds 



(39) 
(40) 



The function / depends on the flux-ratio and the applied 
lens model. 



6.1. Minimizing the negative log-likelihood function 

The user may choose between the three implemented 
methods to minimize the negative log-likelihood function, 
independently of what it looks like. 

1. A simple method to minimize the negative log- 
likelihood function is to calculate its value in a grid, 
specified by start and end values for each parameter 
and the number of grid steps in each dimension. In the 
output file (see Appendix B) the result array is saved 
as an ASCII vector. 

2. The basic idea of the Davidon Variance Algorithm 
(Davidon, 1968) is to calculate the covariance matrix 
by an iterative algorithm. The matrix is obtained by 
using only the function values and the gradient, and 
the minimum value is calculated simultaneously as the 
algorithm converges. This method is much faster in 
finding the minimum than the grid search, the disad- 
vantage however, is that a parabolic shape is assumed 
for the negative log-likelihood function and the esti- 
mated values for the covariance matrix are somewhat 
misleading if this is not the case. 

3. Powell's Quadratically Convergent Method (Powell, 
1964) is a multidimensional minimization algorithm 
that does not require any derivatives. Basically it is a 
method to choose in an efficient way a direction in the 
minimization space and then use linear minimization 
in that direction. This method will however not calcu- 
late any covariance matrix, only the minimum point of 
the function. 



7. Examples of use 

As an example of the usage of SNOC, we show how to 
create a data sample of 2000 Type la SNe with a uniform 
distribution 0.1 < z < 2 and fit the parameters Hm and 
ax (assuming a flat geometry) by performing a maximum 
likelihood analysis. 

In the SNOC input file (see Appendix A), we set: 
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NUMSNE 


2000 




GENDIS 


'uni' 




MAXZ 


2.0 




MINZ 


0.1 




SNTYPE 


'la' 




INTSIG 


-1.0 


(default intrinsic magnitude spread) 


FILTER 


'BZ' 


(redshifted B filter) 


HUBBLE 


0.65 




OMEGAM 


0.3 




DMEGAX 


0.7 




EDSTYPE 


'con' 




ALPHAXO 


-1.0 


(cosmological constant) 


REJTYPE 


3 


(keep only primary images) 


DOFRAC 


1 


(fractional distribution of lenstypes) 


FRACPOI 


0.2 


(20 % point-masses) 


FRACNFW 


0.8 


(80 % NFW halos) 



All other logical options should be set to 'false' or '0'. We 
have used a "standard" cosmology with Q^f = 0.3 and a 
pure cosmological constant fl\ = 0.7 and computed lens- 
ing effects with 20 % point-masses and 80 % NFW halos. 

In Fig. 2, we show the Hubble diagram obtained by 
plotting the various quantities of interest from the SNOC 
output ASCII file (see Appendix B). In the upper panel, 
the luminosity distance, d^, and the lookback time ii.t>. 
(in units of c/Hq) are given as a function of redshift by 
plotting dL and lookb vs zs. 

In the lower panel, the full line is the theoretical 
magnitude-redshift relation for a typical Type la SNe ob- 
served at maximum luminosity in the restframe B band. 
This is obtained by plotting BZmag-intsig vs zs. In the 
same panel the magnitude of each individual SNe is shown, 
including the intrinsic dispersion and gravitational lensing 
effects, i.e., we have plotted BZmag+linsdm vs zs. 

We can now use the SNOC analysis tools to fit, e.g., 
Q, m and ax , assuming Om + ^x = 1 ■ 

In the input for the analysis tool (Appendix A), we 

set: 



2000 
2 

(assuming Qm + ^x = 1) 
-1.0 0.0 200 -1 
0.0 0.4 200 0.3 



nrsne 

nrpar sigma 
flat 
w_0 

omega_M 
NUISM 
nointsig 
gauss 

nodavidonmin 

gridsearch 

nopowellmin 

Fig. 3 shows the confidence contours produced from 
this three-parameter fit, where M has been treated as 
a nuisance parameter. The solid contour that is centred 
around the input values £Im — 0.3, and ax — —1.0 is the 
fitted confidence region when the lensing dispersion is not 
considered, i.e., the keyword nolensef f ect has been spec- 
ified in the input file. The dashed contour shows instead 



- 4 



^ 3 - 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 



26 r 
24 '- 
22 '- 
20 7 
M 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 



Fig. 2. In the upper panel the luminosity distance, d,L, and the 
lookback time ti.h. (in units of c/Ho) are given as a function of 
redshift. In the lower panel the theoretical magnitude-redshift 
relation for Type la SNe (full line) is plotted together with 
the magnitude of each individual SNe including the intrinsic 
dispersion and gravitational lensing effects. 



-0.6 



a* -os 




0.05 0.1 



0.15 0.2 0.25 0.3 0.35 0.4 



Fig. 3. Three-parameter fits in the (S7m , a^-plane using the 
data set shown in Fig. 2. The contours show the 1.51<r level, 
which corresponds to a probability content of ~ 68 % of in- 
cluding both parameter in the region. 

the fitted result when the addlensef f ect parameter has 
been specified in the input file and the lensing effects are 
considered. As can be seen, this will induce a significant 
bias in the parameter determination, an effect discussed 
in further detail in Amanullah et al. (2002). Both contours 
shows the 



min(£) 



1.51 2 
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level of the negative log-likelihood surface, C, defined in 
Eq. (37). This level corresponds approximately to a con- 
fidence region of ~ 68 % of including the true value of the 
two parameters of interest. 

8. Summary 

As we enter a phase of high-precision cosmological mea- 
surements, sophisticated tools for data analysis are re- 
quired, especially to evaluate potential systematic effects 
related to the method. The Monte-Carlo simulation pack- 
age SNOC is mainly designed for the study of the preci- 
sion and possible biases of measured cosmological param- 
eters from high-z Type la SNe. Besides calculating the lu- 
minosity distance for arbitrary cosmological parameters: 
iJ , ft M ,Q x ,c<x(z), the code also allows for the estimates 
of observed magnitudes for high-z SNe taking into account 
the inhomogeneity of matter and dimming by, e.g., dust 
or more exotic process such as photon-axion oscillations. 
Along with the simulation, a maximum-likelihood analy- 
sis package has been developed for fitting of cosmological 
parameters and possible bias. 
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Appendix A: Summary of SNOC input parameters 
The Input for the Monte-Carlo tool 



parameter name description 



NUMSNE 


number of SNe to be generated 


SEED 


seed for random number generator (if negative use computer clock) 


SNTYPE 


SN type 


Redshift distribution 


GENDIS 


functional form of SN redshift distribution option 


MEANZ 


mean redshift of SN sample (for Gaussian distribution only) 


SIGMAZ 


a in redshift of SN sample (for Gaussian distribution only) 


MAXZ 


maximum SN redshift 


MINZ 


minimum SN redshift 


ZFILE 


name of file with tabulated z-distribution of SNe (optional) 


TAUIA 


r-parameter to in SFR redshift distribution (Dahlcn & Fransson, 1999) 


Lightcurve 


INTSIG 


intrinsic magnitude spread of SNe 


BUILDLC 


option to build SN lightcurves 


LCPOINT 


number of points in lightcurve (optional) 


LCNOISE 


measurement error in each lightcurve point 


LCSPACE 


number of days between lightcurve points 


LCTMIN 


earliest day in lightcurve 


FILTER 


default filter for lightcurve (BVRIJ) 


USFILT 


user specified filters 


STRETCH 


lightcurve "stretch", relevant for Type la SNe (Perlmutter et al., 1997) 


STRCORR 


lightcurve "stretch-brightness" correlation (Perlmutter et al., 1997) 


STRSIG 


intrinsic lightcurve "stretch" standard deviation 


Cosmological parameters 


HUBBLE 


Hubble parameter, h (H = /i-100 km s _1 Mpc -1 ) 


OMEGAM 


mass density, Qm(z = 0) 


OMEGAX 


dark energy, fix (z = 0) 


EOSTYPE 


dark energy equation of state form (constant or time varying) 


ALPHAXO 


dark energy equation of state parameter, w = Px / px at z=0 


ALPHAX1 


dark energy equation of state parameter, first derivative, w\ = 


EOSFILE 


dark energy equation of state parameter, tabulated file with w(z) (optional) 


Gravitational lensing 


LENS 


parametrization of matter lenses (compact objects, uniform, SIS or NFW) 


DOFRAC 


option to use alternating lens types 


FRACHOM 


fraction of path in homogeneous matter 


FRACP0I 


fraction of path in an environment of compact objects 


FRACUNI 


fraction of path in an environment of uniform objects 


FRACISO 


fraction of path in an environment of SIS mass distributions 


FRACNFW 


fraction of path in an environment of NFW mass distributions 


CELLSZ 


initial cell size [Mpc] (when MTYPE 1) 


MATTERB 


size of homogeneous matter ball [Mpc] (when LENS 'uni') 


MULTI 


option to treat multiple images of gravitationally lcnscd SNe 


REJTYPE 


event rejection level for secondary (lensed) images 


CLUSTER 


option to simulate lensing due to foreground galaxy cluster 


ZCL 


redshift of foreground galaxy cluster 


VELDISP 


dispersion velocity of foreground galaxy cluster 


MASSCLU 


mass of foreground galaxy cluster 


RADCLU 


radius of foreground cluster image field 


Dust extinction 


DOGREY 


option to add intergalactic dust in the line-of-sight 


MINWAVE 


shortest observable wavelength 
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MAXWAVE 


longest observable wavelength 


NUMPOI 


number of points in the wavelength range 


GREYTAB 


option to use tabulated grey dust parameters vs redshift 


GDFILE 


file with grey dust tabulated data 


RV 


intergalactic dust parameter, see Eq. (23) 


GREYSC 


mean free path for intergalactic dust scattering 


ZLIM 


turning redshift in dust model, see Eq. (25) 


PI 


exponent in dust density redshift evolution below ZLIM, see Eq. (25) 


P2 


exponent in dust density redshift evolution above ZLIM, see Eq. (25) 


DOHOST 


option to simulate host galaxy type, brightness and dust extinction 


RVHOST 


Ry for host galaxy dust 


HOSTSC 


mean free path for host galaxy dust scattering 


DOGLXY 


option to consider dust extinction from foreground galaxies 


RVGLXY 


Ry for foreground galaxy dust 


GLXYSC 


mean free path for foreground galaxy dust scattering 


Photon-axion oscillations 


DOMIX 


option to consider photon-axion oscillations 


BO 


mean comoving magnetic field 


BOWID 


width of uniform dispersion of magnetic field 


MMA 


strength of photon-axion coupling 


MA 


axion mass 


NE 


mean comoving electron density 


NEWID 


width of uniform dispersion of electron density 


K-corrcctions 


DOKCORR 


calculate K-corrections at maximum from restframe B to observed BVRIJ 


KFULL 


write out magnitude distortions from host galaxy dust or line-of-sight dust 


Galaxy parameters 


EVOLUM 


evolution if host galaxy luminosities according to L(z) = L(0)(l + Z ) KVULUM 


MTYPE 


equal mass galaxies or weighted by Schechter function 


SCHEAL 


Schcchtcr function parameter 


SN search 


DOSEAR 


option to simulate search conditions 


TIMEGAP 


time between reference and discovery images in search mode 




Table 1: The SNOC Monte-Carlo tool steering parameters 



The Input for the Analysis Tool 

parameter name description 

nrsne number of SNe to use in the fit 

zrange only use SNe within specified z-range 

nrpar sigma the number of parameters that should simultaneously lie within a 

specific deviation range, 

flat assume that the Universe is flat 

nof lat do not assume that the Universe is flat 

Fit parameters 

These parameters specify which of the cosmological parameters that should be fixed 
and which that should be fitted. Each line consist of a parameter name and four 
numerical values, s, e, g and c. For a parameter to be fitted s and e cannot be 
identical, and g should be an integer greater than 1. For a grid search s and e specifies 
the grid interval for each parameter and g then number of steps. The value c should 
be the correct value of each parameter that was used in the Monte-Carlo. For the 
iterative minimization methods a start value will randomly be selected between s and 
e. 
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w_0 s e 


g 


c 






W_l S 6 


g 


c 






omega_M 


s 


e 


g 


c 


omega_X 


s 


e 


g 


c 


mscript 


s 


e 


g 


c 


hubbleO 


s 


e 


S 


c 



NUISM 
NUIFP 



non-redshift dependent part of the equation of state parameter 
linearly redshift dependent part of the equation of state parameter 
mass density, Q-m (z = 0) 
dark energy, fix (z = 0) 

intercept in the Hubble diagram, A4, only used for Hubble dia- 
gram fits 

Hubble parameter, h (Hq = h- 100 km s" 1 Mpc^ 1 ), only used for 
fits of multiply imaged objects 

Nuisance parameters 
fit Al as a nuisance parameter 
fit the fraction of compact objects in the Universe 



Type la fits 



gausserror 

addintsig 

nointsig 

addlensef f ect 

nolensef f ect 

magresiduals 

nomagresiduals 

gauss 
gausslin 



the uncertainty in the magnitude a 
add a dispersion to the measured magnitudes 
do not add a dispersion to the measured magnitudes 
add the magnitude dispersion due to lensing 
do not add the magnitude dispersion due to lensing 
calculate the magnitude residuals from the best cosmology fit 
do not calculate the magnitude residuals from the best cosmology 
fit 

assume Gaussian magnitude distribution 

assume a modified Gaussian magnitude distribution to handle 
gravitational lensing 

Multiple image fits 

domulti do a multiple image fit instead of the default Type la fit 

ibandmax maximum allowed I-band magnitude for the SNe 

lenstype which lens model to use, SIS or POI 

etimedelay uncertainty in the time delay in days 

eimsep uncertainty in the image separation in arcsecs 

imrat uncertainty in the flux ratio 

ezl uncertainty on lensing galaxy redshift 

emodel systematic error on lensing model 

Minimization method 
davidonmin use the Davidon Variance Algorithm 

nodavidonmin do not use the Davidon Variance Algorithm 
gridsearch use the grid search 

nogridsearch do not use the grid search 
powellmin use Powell's Quadratically Convergent Method 

nopowellmin do not use Powell's Quadratically Convergent Method 

Table 2: SNOC analysis steering parameters 



Appendix B:The SNOC Output 

The Output of the Monte-Carlo Tool 

SNOC produces an ASCII output file that can be easily interpreted and piped into histogramming tools. 6 Table 3 
shows the first event record of a SNOC simulation. 

SNOC output file: example 

variable value I description 



We use the CERNLIB package PAW for that purpose 
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event 


1 


event number 


sntype 


la 


OAT 

bJN type 


nubble 


0.6500 


TT 111 

Hubble parameter 


omegam 


0.2800 


tt M (z = 0) 


omegax 


0.7200 


Mx(z = 0) 


alphaxO 


-1.0000 


w (z = 0) 


Ntotal 


237 


l r pat l 1 1 

number ol blNe sq.deg year 


nzbins 


28 


number of z-bins for rates 


zcent 


0.125 0.175 0.225 0.275 0.325 0.375 


center of z-bin 




0.425 0.475 0.525 0.575 0.625 0.675 






0.725 0.775 0.825 0.875 0.925 0.975 






1.025 1.075 1.125 1.175 1.225 1.275 






1.325 1.375 1.425 1.475 




zrate 


0.749 1.414 2.154 2.770 3.741 4.576 


l * rr 1 • i hit 1 

differential SN rate 




5.076 6.144 6.896 7.161 8.218 8.818 






8.828 9.824 10.260 10.043 10.961 






11 ocn 1 n oca 1 1 /^ooii c/2 1 1 oin 

11.250 10.850 ll.ooo ll.ooo 11.319 






10 AGO 10 1 T£! 11 TOO IO OOO 10 OOO 

12.083 12.156 11.523 12.223 12.223 






11.527 




zs 


0.4318 


redsniit or blN (m event #=1J 


~i 1,"U 

lookb 


U.OZD4 


lookback time (Hubble time units) 


lookGyr 




lookback time in Giga-years 


uL 


U.OOoD 


luminosity distance in units of 


BZmag 


22.9903 


effective B-band magnitude of event 


intsig 


0.2662 


intrinsic B-band offset from mean B-band la magnitude 


emptyb 


1.0241 


luminosity ratio empty beam/filled beam 


emptydm 


0.0259 


*xl iX-i- .1-1 -G11J1 

magnitude offset empty beam - filled beam 


lensrat 


1.0212 


lensing flux ratio demagnification 


linsdm 


0.0228 


lensing demagnification (magnitude) 


ncells 


000 


number of galaxy cells along path 


ncaustics 





number of parity-inversions of image along path 


nre j ect 


n 
U 


number of image rejections before keeping event 


ttot 


1508.5641 


travel time (Mpc) 


galtype 


obc 


host galaxy type 


galabsB 


-10. 4y© / 


host galaxy absolute magnitude 


galBband 


26.2917 


host galaxy B-band magnitude 


gal Vb and 


25.0cio2 


host galaxy V-band magnitude 


galRband 


24.2002 


host galaxy R-band magnitude 


gallband 


23.5455 


host galaxy I-band magnitude 


hostdm 


a aaaa 
U.UUUU 


effective B-band extinction in host 


greysc 


COOAA AAA 

58200.000 


intergalactic dust mean free path (Mpc) 


r_v 


9.5000 


Rv intergalactic dust 


zgreylim 


0.5000 


intergalactic dust, density turn-point, see h,q. (25) 


greyexpl 


3.0000 


intergalactic dust, exponent below turn-point , see hq. (25) 


greyexp2 


0.0000 


intergalactic dust, exponent above turn-point, see tq. (25) 


greydm 


0.5118 


intergalactic dust magnitude extinction ol effective B-band 


numgrey 


30 


number of wavelength bins for intergalactic reddening 


lambdagr 


4000.0 4448.3 4896.6 5344.8 


wavelength central value 




5793.1 6241.4 6689.7 7137.9 7586.2 






0U04.U o4:OZi.o oyoi.u yo ( y.o uozi i .u 






10275.9 10724.1 11172.4 11620.7 






12069.0 12517.2 12965.5 13413.8 






13862.1 14310.3 14758.6 15206.9 






15655.2 16103.4 16551.7 17000.0 
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greyvec 



mixdm 

KcorBBO 

KcorBVO 

KcorBRO 

KcorBIO 

KcorBJO 

KcorBBl 

KcorBVl 

KcorBRl 

KcorBIl 

KcorBJl 

KcorBB2 

KcorBV2 

KcorBR2 

KcorBI2 

KcorBJ2 

KcorBB3 

KcorBV3 

KcorBR3 

KcorBI3 

KcorBJ3 

KcorBB4 

KcorBV4 

KcorBR4 

KcorBI4 

KcorBJ4 

KcorBB 

KcorBV 

KcorBR 

KcorBI 

KcorBJ 

stretch 

strcorr 

npoints 

lcdate 

lcmag 

lcnoise 
dayref 
daynew 
magdif f B 
magdif f V 
magdif f R 
magdif f I 
magdif f J 
magnewB 
magnewV 



0.558E-01 
0.532E-01 
0.476E-01 
0.415E-01 
0.356E-01 
0.300E-01 
0.252E-01 
0.214E-01 
0.185E-01 
0.161E-01 0. 



0.560E-01 
0.514E-01 
0.456E-01 
0.395E-01 
0.337E-01 
0.283E-01 
0.239E-01 
0.204E-01 
0.176E-01 
154E-01 



0.539E-01 

0.548E-01 

0.496E-01 

0.436E-01 

0.375E-01 

0.318E-01 

0.267E-01 

0.226E-01 

0.194E-01 

0.168E-01 

0.2443 

1.2692 

-0.1974 

-0.7003 

-0.7416 

-0.1928 

1.2808 

-0.1873 

-0.6916 

-0.7349 

-0.1892 

1.2692 

-0.1974 

-0.7003 

-0.7416 

-0.1928 

1.3246 

-0.1434 

-0.6500 

-0.6978 

-0.1665 

1.4405 

-0.0273 

-0.4390 

-0.4464 

0.1988 

1.5257 

0.0546 

-0.3646 

-0.3843 

0.2350 

0.9457 

0.0543 



-13.54 -6.77 0.00 6.77 13.54 20.31 
27.08 33.85 

0.42 -0.05 -0.58 -0.28 0.13 0.54 1.15 

1.61 

0.1000 

-12.1837 

3.1770 

-0.9194 

-1.2900 

-1.4338 

-1.5398 

-1.3374 

1.5924 

-0.0228 



differential extinction 



effective B-band extinction due to photon-axion oscillations 
KeB-correction, LC maximum 
Key-correction, LC maximum 
KeR-correction, LC maximum 
Kei-correction, LC maximum 
Kej-correction, LC maximum 
KeB-correction (includes host galaxy extinction) 
Key-correction (includes host galaxy extinction) 
KeR-correction (includes host galaxy extinction) 
Kei-correction (includes host galaxy extinction) 
Kej-correction (includes host galaxy extinction) 
KBB-correction (includes foreground galaxy extinction) 
Key-correction (includes foreground galaxy extinction) 
KeR-correction (includes foreground galaxy extinction) 
Kei-correction (includes foreground galaxy extinction) 
Kej-correction (includes foreground galaxy extinction) 
KeB-correction (includes intergalactic extinction) 
Key-correction (includes intergalactic extinction) 
KeR-correction (includes intergalactic extinction) 
Kei-correction (includes intergalactic extinction) 
Kej-correction (includes intergalactic extinction) 
KeB-correction (includes photon-axion oscillation) 
Key-correction (includes photon-axion oscillation) 
KeR-correction (includes photon-axion oscillation) 
Kei-correction (includes photon-axion oscillation) 
Kej-correction (includes photon-axion oscillation) 
KeB-correction total extinction folded 
Key-correction total extinction folded 
KeR-correction total extinction folded 
KBi-correction total extinction folded 
KBj-correction total extinction folded 
lightcurve stretch 
stretch-brightness correction 
number of lightcurve points 
dates from maximum 

nifl(i) - BZmag (R-band lightcurve) 

applied noise standard deviation in LC points 
days from maximum for reference image (rest frame) 
days from maximum for search image (rest frame) 
SN B-band magnitude difference new - ref epoch 
SN V-band magnitude difference new - ref epoch 
SN R-band magnitude difference new - ref epoch 
SN I-band magnitude difference new - ref epoch 
SN J-band magnitude difference new - ref epoch 
SN B-band magnitude (-BZmag) at discovery epoch 
SN V-band magnitude (-BZmag) at discovery epoch 
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magnewR -0.6042 SN R-band magnitude (-BZmag) at discovery epoch 

magnewl -0.7418 SN I-band magnitude (-BZmag) at discovery epoch 

magnewJ -0.0819 SN J-band magnitude (-BZmag) at discovery epoch 

end end of event # 1 

Table 3: Extract of SNOC output showing the first event record 



The Output of the Analysis Tool 

SNOC output file: example 

# THIS IS THE OUTPUT GENERATED BY SNALYS 

# ====================================== 

# The number of supernovae used: 2300 

# The number of parameters estimated: 2 

# Flat universe assumption: no 

# Execution started at: Wed Mar 6 16:24:15 2002 

# Total execution time: 00:01:33.00 

# Used CPU-time: total = 00:01:35.43 

# user = 00:01:35.34 

# system = 00:00:00.07 



# Powell Minimization 

# =================== 

# Start values 
omega_m . 333 
omega_x 0.733 

# The min of the ML function (powell) : 0.010 

# Cosmology (powell) 
omega_m 0.299860 
omega_x 0.699683 

# Davidon Minimization 

# ==================== 

# Start values 
omega_m . 300 
omega_x 0.700 

# The min of the ML function (davidon): 0.010 

# Cosmology (davidon) 

omega_m 0.299983 +/- 0.023836 
omega_x 0.699974 +/- 0.046131 

# Correlation matrix (davidon) 

omega_m omega_x 
omega_m 1.00000 0.95052 

omega_x 0.95052 1.00000 



# Grid Search 

# =========== 

# The min of the ML function (grid): 6.729 

# Cosmology (grid) 

omega_m 0.333333 + 0.33E-01 - 0.67E-01 
omega_x 0.733333 + 0.00E+00 - 0.67E-01 

# Conditions for grid search. 

# Parameter start stopnr.step 

omega_m 0. 10 0.40 10 
omega_x 0.40 1.00 10 

# The ML-function from Grid Search 



16 



Goobar et al.: SNOC: a Monte-Carlo simulation package for high-z supernova observations 



#omega_m omega_x 
82.4738235 
42.1675644 
51.1707954 



